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Abstract 

In this work we study the O + (1540) in the framework of QCD sum 
rules based on (ud) 2 s diquark clustering as suggested by Jaffe and 
Wilczek. Within errors, the mass of the pentaquark is compatible 
with the experimentally measured value. The mass difference between 
the © + and the pentaquark with the quantum numbers of the nucleon 
amounts to 70 MeV, consistent with the interpretation of the iV(1440) 
as a pentaquark. 
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Recently, several experiments [1-10] have observed a new baryon resonance 
+ (154O) with positive strangeness. Therefore it requires an s and has a 
minimal quark content of five quarks. This discovery has triggered an in- 
tense experimental and theoretical activity to clarify the quantum numbers 
and to understand the structure of the pentaquark state. The O has the 
third component of isospin zero and the absence of isospin partners sug- 
gests strongly that the is an isosinglet what we also assume in this work. 
A puzzling characteristics of the 6 is its narrow width below 15 MeV. A 
suggestive way to explain the small width is by the assumption of diquark 
clustering. The formation of diquarks presents an important concept and 
has direct phenomenological impact [11]. Two models have been proposed 
based on the strong attraction of the (-U(i)-diquarks: one by Karliner and 
Lipkin [12] where the pentaquark is described as diquark-triquark system 
in a non-standard colour representation. The other one is due to Jaffe and 
Wilczek [13, 14] and describes the G as bound state of an s with two highly 
correlated (wci)-diquarks. In this work we investigate the second approach by 
Jaffe and Wilczek in the framework of QCD sum rules. In principle, as was 
discussed in [15], even a mixing between the two states could be possible. 
However, an estimation of such a potential mixing would require a detailed 
investigation of the model by Karliner and Lipkin. The basis of the sum rules 
was laid in [16] and their extension to baryons was developed in [17]. The 
assumptions of the model are incorporated by an appropriate current. Since 
the sum rules are directly based on QCD and keep the analytic dependence 
on the input parameters, they can help to differentiate between the models 
and to test their features. The relevance of the diquark picture within the 
context of the sum rules was shown in [18]. Several sum rule investigations 
for the pentaquark already exist [19-22] which, however, are based on dif- 
ferent models or currents. The diquark models for the pentaquark have also 
been investigated within other approaches [23]. 

In the model by Jaffe and Wilczek the (w<i)-diquarks have zero spin and 
are in a 3 C and 3/ representation of colour and flavour. In order to combine 
with the antiquark into a colour singlet, the two diquarks must combine 
into a colour 3. The diquark-diquark wavefunction is antisymmetric and 
has angular momentum one. This combines with the spin of the s to total 
angular momentum 1/2 and results in positive parity. In [13] it was suggested 
to interpret the Roper resonance ^(1440) as {ud) 2 d pentaquark state and we 
will study this resonance at the end of our analysis. 

The basic object in our sum rule analysis is the two-point correlation 
function 
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where r)(x) represents the interpolating field of the pentaquark under inves- 
tigation. 

The diquarks have a particularly strong attraction in the flavour anti- 
symmetric J p = + channel. Thus the current contains two diquarks of the 
form 

Q>(x) = e abc Q ab (x) = e abc [u T a C l5 d b ] (x) . (2) 

C denotes the charge conjugation matrix. The two diquarks must be in a 
p-wave to satisfy Bose statistics. Therefore the current contains a derivative 
to generate one unit of angular momentum. The diquarks couple to a 3 C in 
colour to form the current 

v ( x ) = (e abd 5 ce - e abc 5 de ) [Q ab (D»Q cd ) - (D»Q ab ) Q cd ] ^,Cs T e , (3) 

where the covariant derivative for the 3 C is given by D» = d»-ig\\A» l [14]. 
The parity is positive. This current has a different structure than the cur- 
rent in [21] which contains no derivative to produce the angular momentum 
between the diquarks. Inserting the current and neglecting higher orders in 
the strong coupling constant the correlator is given by 

U(x) = (0|T{t7(x)77(0)}|0> = h^CS£ T (-x)C7y 7B l t£(x) , 



Tf£(x) = (e abd 8 ce - e abc S de ) (e a ' b ' d 'S c ' e ' - e a ' b ' c '5 d ' e ' 



x { (i 5 S w (x) l5 CSj a , (x)C) { lb S dd> {x) lb CS T cc , (x)C) 

+ (75^ {x) lb CSl, (x)C) ( lb S dV {x^CSl, (x)C) 

- (ibSm {xhsCSl {x)C lb S dV (x) l5 CSj a , (x)C) 

- (^S bb ix) l5 CSjAx)C l5 S dd ix) l5 CSUx)C)}, (4) 

with = jdx^ + 

q(P) jQ x ii an( j U pp er co l our index indicates the 

propagator on which the derivative is acting. S(x) and S^ s \x) are the light 
and strange quark propagators, respectively. The quark propagator has been 
evaluated in the presence of quark and gluon condensates in [21,24,25], 
where the explicit expressions can be found. Using the following Lorentz 
decomposition for Tf£' = 5 ee 'T^ u /3, 

T»v = 9^fi(x 2 ) + x il x v f 2 {x 2 ) , (5) 

the functions fi(x 2 ) and f2(x 2 ) are determined to 

2 576 240m;; 2Am 4 q - M<K 2 m q (qq) + ^f(^FF) 

) ~ ^12 + ^10 
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12m 3 q (qq) - Am q (g s qaFq) - I6rr 2 (qq) 2 / 6 



7T U X C 



. 1152 576m 2 -48m* + 2567r 2 m„ (gg) - ^f(^FF) 



7T 8 X 16 7T 8 X 14 7r 8 a; 12 
-32m;j(gg) + 32m q (g s qaFq) + 1287r 2 (gg) 2 

7T 6 S 10 



+ 0(l/x 8 ). (6) 



The colour non-diagonal part of Tj^ vanishes for the considered orders. In 
momentum space the correlator can be parametrised as 

n(p) = yn^( P 2 ) + n (1 V)- (7) 

To obtain the phenomenological side we insert intermediate baryon states 
with the corresponding quantum numbers. The matrix element of the is 
parametrised by 

(0|77(0)|e(p)) = / e -«(p). (8) 

Since no experimental information on higher pentaquark states is available 
we make the assumption of quark-hadron duality and approximate the higher 
states by the perturbative spectral density above a threshold s . In fact, the 
uncertainty on so will be one of the dominant errors in the sum rule analysis. 

In order to suppress the higher dimensional condensates and to reduce the 
influence of the higher resonances we employ a Borel transformation defined 
by 

B M = lim t^f-fLV, M 2 = ^ fixed, (9) 
Q2,n^oo r(n) \dQ 2 J n 

with Q 2 = —p 2 . As in [20, 21] we now concentrate on the chirality even 
part in eq. (J2J) which contains the leading order term from the operator 
product expansion. The spectral density p(s) = l/7rlmll^(s + ie) has the 
form 

p(s) = a 6 s 6 + a 5 s 5 + a 4 s 4 + a 3 s 3 + . . . (10) 

The coefficients a« can easily be obtained from the results of eqs. (jU and © 
by inserting the strange quark propagator and performing a Fourier trans- 
formation. The theoretical moments are then given by 

fl(M 2 ) = E M II(Q 2 ) = J^ds P -^e- s ' M2 = a 6 r(7)(M 2 ) 6 

+a 5 r(6)(M 2 ) 5 + a 4 r(5)(M 2 ) 4 + a 3 r(4)(M 2 ) 3 + . . . (11) 

Transferring the continuum contribution to the theoretical side and taking a 
logarithmic derivative with respect to — 1/M 2 , one obtains the sum rule for 
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the mass of the pentaquark, 

k Y!a 6 - k T(8-k)(M 2 ) s - k E 7 _ k 

< = S > ( 12 ) 

£ a 6 _ k T(7-k)(M*y-KE^ k 

k=0 

where E a = 1 - T(a + 1, s /M 2 )/r(a + 1). 

A basic input for the sum rule analysis is the Borel parameter M. The 
sum rule should be stable with respect to M to allow a reliable determination 
of the pentaquark mass. For large values of M the operator product expan- 
sion converges well, however, for small M the expansion becomes problematic 
and thus we restrict the range of the Borel parameter to M > 1.6 GeV. 
Small M suppress the phenomenological continuum part which becomes 
very dominant for large M. Therefore we employ a sum rule window of 
2.5 GeV 2 < M 2 < 4.0 GeV 2 . 

As input parameters in our analysis we use m s = 0.15 GeV, (qq) = 
-(0.267 ± 0.018 GeV) 3 , (ss) = (0.8 ± 0.2) (qq), (g s saFs) = M 2 (ss) with 
Mq = (0.8 ± 0.2) GeV 2 , and (^FF) = 0.024 ± 0.012 GeV 4 [26]. For the 
continuum threshold we use a central value of so = (1.54 + 0.35 GeV) 2 . 
Thus the continuum starts 350 MeV above the measured pentaquark mass. 
This difference should roughly correspond to one radial excitation [20] and 
represents a typical value for sum rule analyses with light quarks as degrees 
of freedom [16]. Fig. ^ shows the mass as a function of the Borel parameter 
M 2 . The sum rule has a good stability with respect to M. As central 
value for the pentaquark mass we obtain me = 1.64 GeV. The two most 
important sources of the error are the choice of the continuum threshold and 
the convergence of the operator product expansion. Since we have substituted 
the phenomenological spectral density, using the assumption of quark hadron 
duality, by the perturbative expansion, the uncertainty on s reflects the 
missing knowledge of the experimental cross section for higher energies. To 
estimate the error on me we vary so between 3.0 GeV 2 < so < 4.1 GeV 2 . In 
fig. [T]we have also plotted the change of with the continuum threshold 
from which we obtain an error of Aitlq w 125 MeV. More phenomenological 
information would be essential to reduce this kind of error. To estimate the 
dependence of the sum rules on the OPE we successively remove the different 
orders. Fig. El shows the convergence of the pentaquark mass including the 
condensate contributions up to a specific power. The inclusion of the higher 
condensates lowers the mass. Using only the leading order perturbative result 
the central value is about 100 MeV larger than the full result. We have 
not included an extra graph for the term cx as since this contribution is 
proportional to the light quark masses and their influence on the analysis 
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Figure 1: m@ as a function of the Borel parameter M 2 for different Sq 
3.5 GeV 2 (solid), s = 4.1 GeV 2 (dotted) and s = 3.0 GeV 2 (dashed). 
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Figure 2: m e for different orders of the OPE, using only the leading order 
perturbative expansion (dotted), with dimension 4 condensates (dashed) and 
including the condensates of dimension 6 (solid). 
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Figure 3: Mass difference m© — tun for different values of the continuum 
threshold, the solid, dashed and dotted lines are for s e — 3.5 GeV 2 and 
s 0N = 3.2 GeV 2 , s © = 4.1 GeV 2 and s 0N = 3.8 GeV 2 and s © = 3.0 GeV 2 
and s n = 2.7 GeV 2 , respectively. 



can be neglected. The four-dimensional condensates lower the leading order 
result by about 50 MeV and the condensates of dimension 6 by another 50 
MeV. We assume that a reasonable error estimate from the OPE would be 
Am© « 75 MeV. Furthermore, contributions to the error also arise from the 
other input parameters which we vary in the ranges presented above. As it 
turns out, their influence on the value of m© is small compared to the errors 
from the continuum threshold and the convergence of the OPE. Adding the 
errors quadratically our final result reads 



In [13] Jaffe and Wilczek suggested to interpret the Roper resonance as {ud) 2 d 
pentaquark state. One can then perform a similar analysis for the iV(1440) 
as has been done for the O by substituting the s antiquark by a d antiquark. 
As central value for the continuum threshold we choose, as in the + case, a 
value of 350 MeV above the ground state mass. For the error range we use 
2.7 GeV 2 < son < 3.8 GeV 2 . Performing a sum rule analysis for the iV with 
the above given parameters, we obtain a mass of mjv = 1.57 ± 0.15 GeV. 
Similar as it has been done in [21], in fig. El we plot the mass difference 
m© — irtN for different values of the continuum thresholds. The mass split- 
ting between the pentaquark states comes out to be about 70 MeV. The error 



m© = 1.64 ±0.15 GeV. 



(13) 
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represented in fig. |3] is based on the assumption that the continuum thresh- 
olds have the same offset for both pentaquark states. Phenomenologically, 
these values can be different and one should add to the error a part of the 
uncertainty from so given in fig. ^ Thus the error can easily amount to 50 
MeV. Though the mass difference is consistent with the interpretation of the 
iV(1440) as a pentaquark, the uncertainty remains large and a reduction of 
the error would be essential to clarify the situation. 

Recently, in [27] it has been argued that one should subtract all possible 
colour- singlet meson-baryon contributions from the pentaquark current. We 
believe that this claim is not correct. Nothing is wrong to use the current 
of eq. This current contains also 2-particle intermediate states which 

have to be added to the phenomenological side. However, at energies around 
the pentaquark mass we expect the pentaquark contribution to dominate the 
spectral density. Apart from KN production whose threshold lies somewhat 
below the pentaquark energy other intermediate states start at higher energy. 
Therefore it is expected that the baryon-meson continuum contribution only 
becomes important at energies much above the pentaquark mass. In this 
energy range the spectral density is suppressed by the exponential in eq. (|TT|) 
and the correlator should be well approximated by the assumption of quark- 
hadron duality. Furthermore, the current rj(x) is based on the assumption 
of diquark formation. Subtracting partial contributions from the OPE side 
changes the pentaquark current and can remove contributions relevant for 
the diquark formation. Thus these contributions can form an important part 
of the pentaquark and should not be subtracted. 

To summarise, we have performed a QCD analysis based on the approach 
by Jaffe and Wilczek. We obtain a sum rule that is stable over the Borel 
parameter M and reproduces the mass of the pentaquark within errors. The 
error is to a large part due to the lack of experimental information above the 
pentaquark energy. Furthermore, a complete calculation at next-to-leading 
order would help to quantify the uncertainties in the theoretical expansion. 
However, with the complex structure of the current and given the fact that 
this includes a calculation of five loops, this is a difficult task. We have 
also performed an analysis for the pentaquark with the quantum numbers of 
the nucleon and have shown that the interpretation of the Roper resonance 
iV(1440) as (ud) 2 d pentaquark state is consistent with the sum rules. It is 
important to note that the sum rules are directly based on QCD and thus, 
apart from the structure of the current, do not contain further model as- 
sumptions. It would be interesting to see if lattice calculations could confirm 
these findings. First lattice calculations exist [28] which, however, are based 
on different interpolating currents and whose results are not yet conclusive. 
Further advance in two directions seems feasible: higher lying pentaquark 
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states with different quantum numbers and internal structure could be in- 
vestigated and a QCD analysis based on the approach by Karliner and Lipkin 
should be done. This might help to understand the specific features of the 
models and to differentiate between the approaches. 
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